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Abstract The spatial-temporal distribution of absorp- 
tion-line systems (ALSs) observed in QSO spectra 
within the cosmological redshift interval z =0.0-4.3 is 
investigated on the base of our updated catalog of ab- 
sorption systems. We consider so called metallic sys- 
tems including basically lines of heavy elements. The 
sample of the data displays regular variations (with am- 
plitudes ~ 15 - 20%) in the z-distribution of ALSs as 
well as in the 77-distribution, where 77 is a dimensionlcss 
line-of-sight comoving distance, relatively to smoother 
dependences. The 77-distribution reveals the periodic- 
ity with period A77 = 0.036 ± 0.002, which corresponds 
to a spatial characteristic scale (108 ± 6) h~ x Mpc or 
(alternatively) a temporal interval (350 ± 20) h^ 1 Myr 
for the ACDM cosmological model. We discuss a possi- 
bility of a spatial interpretation of the results treating 
the pattern obtained as a trace of an order imprinted 
on the galaxy clustering in the early Universe. 



Keywords galaxies: quasars: absorption lines ■ 
ies: high-redshifts 



1 Introduction 
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In our previous pape rs (e.g. 

Ryabinkov et al.l 120071 hereafter Paper I) we studied 
the spatial-temporal distribution of absorption-line sys- 
tems (ALSs) imprinted in spectra of quasars (QSOs). 
We regard that the redshifts z a of absorption lines are 
cosmological, and detected ALSs are associated with 
ionized gas in intervening galaxies or clusters of galaxies 
at cosmological distances along lines of sight to origi- 
nal QSOs. Thus a representative sample of ALSs is 
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considered to trace the distribution of matter between 
the observer and QSOs, as well as physical processes 
occurred in different epochs of the cosmological evolu- 
tion. 

In Paper I we explored 2003 ALSs in the redshift 
range z = 0.0-3.7 from our c atalog of absorption sys- 
tems (JRvabinkov et al.l 120031) and restricted ourselves 
by consideration of so called metallic systems includ- 
ing lines of heavy elements. ALSs registered in a spec- 
trum of a quasar might comprise up to 20 - 30 absorp- 
tion lines predominantly within an interval ~ 3000 - 
8000 A. It was shown that the z-distribution of ALSs 
displays a pattern of alternating maxima (peaks) and 
minima (dips) which are statistically significant against 
a smoother dependence (trend). The positions of the 
peaks and dips obtained for different hemispheres on 
the sky turned out to be independent (within statisti- 
cal uncertainties) of observation directions. Moreover, 
it was shown that the sequence of peaks and dips re- 
veals a certain regularity. The power spectrum calcu- 
lated with using a rescaling function t(z) (Eqs. (11) 
and (12) of Paper I) displays the peak at the signifi- 
cance level exceeding 4a relatively to the hypothesis of 
the uniform distribution of ALSs along the axis t(z). 

In this paper we present the results of slightly ex- 
tended statistical analysis performed by similar meth- 
ods as we used in Paper I but in a new bearing. The 
principle point of this paper is the analysis of the ALSs 
distribution in the comoving coordinate system (CS). 
Such a representation enables us to find the periodic- 
ity of the distribution and treat it as a trace of partly 
ordered spatial structures. 

The upper panel of Fig. [T] demonstrates three z- 
distributions Nas(z) of ALSs containing absorption 
lines of heavy elements within the redshift interval 
z = 0.0 - 4.3. One of th em is obtained using 847 sys- 
tems from the catalog bv lJunkkarinen et al.l (|l99ll ). the 
second one is based on 2003 systems registered in the 





Fig. 1 Upper panel: Sliding-average z-distributions of absorption-line systems (ALSs) observed in QSO spectra within the redshift 
interval z = 0.0 - 4.30 with an averaging bin A z = 0.071 and a step of its consecutive shift <5 Z = 0.01; a velocity interval of ALS 
associations (along each li ne of sight) into a s i ngle s ystem Zi (see text) is chosen as Sv = 500km s — x ; thin soli d lines - 847 syst e ms ar e 
taken from the catalog bv ljunkkarinen et al.l 11199 ll ). thin dashed lines - 2003 systems - from the catalog bv lRvabinkov et al.l (2003), 
thick lines - 2322 systems - from the same but updated catalog. Lower panel: same as in the upper panel but for corresponding 



^-distributions (see text) within the appropriate interval r\ 
S„ = 0.005. 



0.0 - 1.87 with an averaging bin A,, = 0.018 and a step of shifting 



spectra of 661 QSO s (emission redshi f ts z e = 0.29 - 4.9) 
from the catalog bv lRvabinkov et al.l ( 20031) . The third 
one is the z-distribution of 2322 ALSs registered in the 
spectra of 730 QSOs (z e =0.035 - 5.01) from the same 
but updated catalog. All three samples of ALSs repre- 
sent the same types of metallic systems including sin- 
gle lines or resonance doublets of ions: C II-IV, Mg II, 
Si II-IV, N V, Al II-III, Fe II etc. They exclude ALSs 
consisting only of neutral hydrogen lines as well as con- 
taining dampe d Lyq absorption system s (DLA). 

Following to iKaminker et al.l ( 20001 ) we use the pro- 
cedure of smoothing away an initial sample of ALS red- 
shifts to reveal possible large scale variations and sup- 
press small scale clustering of the ALSs. All redshifts 
Zj registered in a spectrum of a certain QSO and fallen 
into the velocity interval Sv =500 km s _1 are treated 



as a single absorption system with an averaged redshift 
z a = (Si=i z j)/ n a, where n a is a number of redshifts 
included in the group. Three distributions in the up- 
per panel of Fig. Q] are obtained using so-called sliding- 
average approach which represents a set of consecutive 
displacements of the averaging bin A z = 0.071 along 
z-axis with a step S z = 0.01. 

The lower panel of Fig. [1] represents the same three 
distributions but calculated relatively to a di mension 
less l i ne-of-site comoving distance tj (z) (e 
19931 iKavser et al.lll997riHogglll99 



Harrison 



■n(zi) 



l 



o ^n m (i + z) 3 + n A 



dz, 



(1) 



where i is a numeration of all ALSs i = 1,2, ...N to t sam- 
pled as described above, N to t = 847, 2003, and 2322, 



respectively. We use the ACDM-cosmological model 
with the dimensionless density parameters Q m = 0.23 
and Qa = 1 - ^m = 0.77. 

A comparison of the z-distributions reveals simi- 
lar patterns of the peaks and dips with rather small 
(~ 15 — 20%) amplitude of variations relatively to 
smoother dependences. The positions of the major- 
ity of peaks and dips remain the same after the exten- 
sion of statistics and some of them become more signif- 
icant (see Paper I). We have similar situation with r\- 
distributions, although there is noticeable lack of statis- 
tics at r\ < 0.6. 

The main accent of the present paper is placed on a 
search for a periodicity of the ALSs radial distribution 
with respect to the variable r\. In Section 2 we focus 
on the periodicity of the ^-distribution. In order to 
avoid sensitivity of our results to the choice of an aver- 
aging bin and to the procedure of a trend elimination, 
which are mutually dependent, we employ in Sections [2] 
and [4] a point-like statistical technique. In Section 3 we 
examine the periodicity of the one-dimensional correla- 
tion function calculated in the comoving CS. In Section 
4 we estimate a trace of the same periodicity in the dis- 
tribution of resonance absorption doublets Mg II based 
on the data of the Sloan Digital Sky Survey (SDSS). 
Conclusions and discussion of the results in application 
to the Large-Scale Structure are represented in Section 
5. Two toy models for partly ordered simple-cubic lat- 
tice displaying similar periodicities are discussed in Ap- 
pendix. 



2 Periodicity of ^-distribution 

In this Section we employ so called point-like statisti- 
cal approach at which a sequence of points r\i = rj(Zi) 
calculated with use of Eq. (TTJ is analysed. To sam- 
ple ALSs for this analysis we use hereafter the same 
averaging velocity interval Sv — 500km s _1 for the as- 
sociation of ALSs into a single redshift Zi as in Fig. [TJ 
Additionally, we exclude all ALSs with redshifts z a be- 
longing to a region associated with their host QSO, i.e., 
at \z e — Zo|/(l + z e) < Av/c, where z e is the emission 
redshift of QSOs and Av is the minimal velocity shift 
(along a line of sight) adopted for the sample of ALSs. 
In this paper we choose Av = 1100 km s _1 , thus the 
actual sample for our analysis consists of N to t = 2167 
systems. 

To verify the periodicity of the ^-distribution we 
calculate a power spectrum for the whole sequence of 



points rji 



P(k) = 




(2) 



where Iq > 1 corresponds to a lowest point of the se- 
quence rj(zi), f\ = ^(^Ntot) — Vh i s the whole interval 
under consideration, k is a harmonic number. A peri- 
odicity yields a peak in the power spectrum V = max 
(P(k)) with a confidence probability 



/3=[l-exp(-P)], 



(3) 



which is defined with respect to the hypothesis of the 
Poisson distribution of t]i . 

Following Eq. (J2J we performed additionally numer- 
ous calculations of the power spectra P(k) for dif- 
ferent realizations of randomly distributed points Xi 
(1 < i < N to t) within intervals Da;. We assumed that 
density of points Xi obey to Gaussian distribution with 
equal values of an expectation and a variance deter- 
mined by a trend (smooth dependence). Using Kol- 
mogorov criterion we obtained appropriate agreement 
between distributions of peak values V at various k and 
the exponential distribution. These simulations confirm 
that a significance probability at certain critical levels 
V CI may be estimated as 



P{P(k) > P CI } = exp{-P cr }, 



(4) 



and Eq. §3§ does for estimations. Our simulations have 
shown that the point-like analysis may be treated as 
the limit case of the standard approach with succes- 
sively narrowed independent bins. The presence of the 
trend is a source of most powerful long- waved harmon- 
ics (smallest k) in the power spectra which should be 
rejected in our consideration. 

The power-spectrum analysis of the ALSs is carried 
out with using Eq. [2] for successively reduced intervals 
£min — z — 4-3 (77 m i n < r\ < 1.87), the low boundary 
being shifted upward from z m i n = to z m j n = 2.0 with 
a step Sz m i n = 0.01. Filled regions in the upper panel of 
Fig. [2] indicate positions of the main peaks with signif- 
icance > 3<7 on the (z m ; n ,k)-plane at integer harmonic 
numbers. All marked peaks correspond to the same pe- 
riod A?y = (?7 max - Jy™„)/k = 0.036 ± 0.002. The lower 
panel displays the value of peaks P max (k) calculated for 
all numbers k corresponding to the filled regions on the 
upper panel. Note the presence of periodicity in the 
ry-distribution of ALSs for the whole interval of z = 0.0 
- 4.3 (77 = 0.0 - 1.87) at the significance level > 3a. 
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Fig. 2 The set of peaks of the power spectra Pmax(k) calculated according to Eq. J2J) in a wide region of harmonic numbers k 
at the significance > 3<r; intervals of z (or rj) are limited from below by z m i n (»?min) shifted successively upward along an interval 
0.0 < .z m in < 2.0 with a step <5z m i n = 0.01. The whole sample contains 2167 values m (see text). All visible peaks correspond to the 
same period An = 0.036 ± 0.002. Upper panel indicates positions of the main peaks on the (z m i n , k)-plane; dimensions of filled regions 
qualitatively denote the values (amplitudes) of the peaks, smallest regions (points) symbolize 3<r significance level. Lower panel plots 
the values of the highest-peaks P max (k) versus the integer number k corresponding to the regions marked on the upper panel. The 
horizontal dash lines specify the significance levels 3<r, 4<r, and 5<r calculated with use of Eq. J3]l 



The highest peaks (~ 4.5c) take place for reduced in- 
tervals z min > 1.64 (7j m i n > 1.148) where statistics is 
essentially better. 

For illustration, the upper panels in Figs. [3] and |4] 
demonstrate the power spectra P(k) calculated in a 
wide region of the harmonic number k for the whole 
interval z = 0.0 - 4.30 (77 = 0.0 - 1.87) and for reduced 
one z = 1.64 - 4.30 (77 = 1.148 - 1.87), respectively. 
The lower panel in Fig. [3] displays a relatively sharp 
dependence of the main peak amplitudes P(k=51) on 
the density parameter fl m , which could prevent one to 
bring out the period A77 in the whole interval of 77, e.g., 
at f2 m ~ 0.3 (used in Paper I). We assume that this 
sharp dependence is bounded with a lack of statistics 
a t V 5; 0-6, although it should be verified in further 
analysis. 



The peak at k=20 (A77 = 0.036) in the upper panel 
of Fig. |4] exceeds the significance level 4.5c. The lower 
panel demonstrates the dependence of the peak ampli- 
tude on Q m . One can see that in contrast with the 
lower panel in Fig. [3] the peak values become almost 
independent of f2 m at Sl m > 0.2. 



3 One-dimensional correlation function 

For verification of the periodicity obtained in Sec- 
tion [2] we calculate additionally a two-point corre- 
lation function £,(Srj) for the sample of ALSs. In 
our consideration a variable Srj substitutes the comov- 
ing distance r between pairs of sampled objects in 
the standard two-point correlation function £(r) (e.g., 
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Fig. 3 Upper panel represents the power spectra P(k) calcu- 
lated for the whole interval rj = 0.0 - 1.87 (z = 0.0 - 4.30) at 
f2m = 0.23; the main peak corresponds to k=51 and Ar; = 0.037. 
Lower panel plots the dependence of the peak value P(k=51) on 
l Cm. 



Fig. 4 Same as in Fig. [3] but for the most significant peak 
(4.5 a) in Fig. [2] Upper panel represents the power spectrum 
P(k) calculated for the interval r, = 1.148 - 1.87 (z = 1.64 - 
4.30) at O m = 0.23; the sample contains N A s = 1148 ALSs. 
Lower panel plots the dependence of the main peak P(k=20) on 
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S Mo et all 
13, iHamiltonl 



1992, 



1993, 



Kerscher et al.l 120001 and references therein) . Unlike 
the standard approach the function £(8rj) is based on 
the choice of the single reference center 77 = (z = 0) 
associated with the observer. Accordingly, all radial 
(line-of-site) points r]i(zi) within a layer (j) ± A^/2), 
where A,, is a width of bin, are treated as equivalent 
ones despite of various spatial distances between them. 
In this way we count up all pairs with fixed relative 
line-of-site comoving distance Srjij = \rji — rjj\ between 
arbitrary ALSs numerated by i and j, which are regis- 
tered in various directions. Thus one can write: 



A s i m (dr/) 



(5) 



where A/" bs (5t)) is a number of observed pairs of ALSs 
separated by 5i] = Srjij belonging to the interval 
6rj ± A^/2. The bin width is chosen as A v = 0.018, 
i.e., the same as in low panel of Fig. [TJ Af s im(5v) is 
a number of crossing pairs between the real sample of 
ALSs and the points of a random (Poisson) sample sim- 



ulated in the same interval of 77 with the same smoothed 
distribution (trend) as the real sample. The both sam- 
ples (real and simulated ones) have the same number 
of redshift points. 

Note also that the correlation function ([5]) differs 
from the two-point multicentric correlation function 
which has been calculated for bot h radial (line-of- 
sight) and transverse directions by iGaztahaga et al. 
(I2009allbr) for the spati al distribution of galaxies (see 



also lKazin et al.ll2010bl ). On the other hand, the func- 



tion ([5]) differs from the correlation functions calculated 
by co unting pairs of objects along separate lines of sigh t 
(e.g jQuashnock et alill996l . lBroadhurst fe Jaffe!l2000!) . 
although the results of both approaches seem to be com- 
patible. 

The upper panel of Fig. [5] displays the one-dimensio- 
nal (two-point) correlation function ([5]) calculated in 
units of the appro priate Poisson error cr(£) = (1 + 
£)/\/Xim(<fy) (e.g- IPeacock fc Nicholsor]|l99lh : 



€(6ri) Afobs(Sri) - Kim{Sv) 



*(0 



y/Kim(Sri) 



(6) 
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Fig. 5 Upper panel: two-points correlation function §(<5??)/c(§) (see text) versus the interval of 5n (0.01 < Sn < 0.408) between 
the components of pairs constituted from the same sample of 1148 ALSs as in Fig. [4] Lower panel: power spectrum P/j(k) calculated 
for /i = ^(tj) = £( < 5'?)/ CT (i;) according to Eq. J7J|, k is the harmonic number. 



An interval of values 6r] in Fig. [S] (0.01 < 5rj < 
0.408) is chosen about twice shorter than the 77-interval 
(1.148 < rj < 1.870) to provide its uniform filling by 
various pairs of ALSs. One can see the sequence of 
positive and negative peaks with the significance > 3ct 
with respect to zero level. Let us notice the presence of 
a long-range order in the dependence of £,{5r])/a(!;) on 
J77. We discuss this effect in Section [5) 

The lower panel in Fig. [5] represents the power spec- 
trum calculated for the value /j, = n(Srf) — £(6r))/a(£) 



according to the equation: 



P/xGO 



A^[/(k) 



N u 



2J Mm sin 

m— 1 



^ Mm COS 



27rk Srj.„ 

6fj 



27rk 8i]„ 
6f] 



(7) 



where (x m = 
of points m 



^(Si] m ), the values 5i] m run over a set 
= 1,2, ..JVp of the variable Srj from 0.01 
to 0.408. Here we employ the sliding bins A,, = 0.018 
with centers at Srj m which transit from one center to 
another one by successive shifts h„ = 0.001; A M = 399 is 
the number of the points 5r) m within the whole interval 



Sf) = 0.398. The value U(k) in the denominator of (O 
may be presented as 



f/(k) 



where 



J7(k) = 



U(k) at U(k) > 1 
1 at C/(k) < 1, 



Sm( — j-r- 1 ) 



()f/ 



<5r) 7rkA r/ 



A, 



(8) 



(9) 



It is known from spectral analysis (e.g., lBathlll974l ). 
that the sliding-average procedure corresponds to a low- 
frequency filtration with a rectangular function as a fil- 
ter. Fourier transform of this function is represented by 
sin(w/2)/(cj/2), where uj is a dimensionless Fourier fre- 
quency. Note that in the considered case a factor A v /h v 
in Eq. © is equal to 18. Our special simulations of 
random catalogs of points have shown that the normal- 
ization factor (|HJ| , ([9]) in the denominator of expressions 
analogous to Eq. (|7|) compensates a low-frequency dis- 
tortion through the sliding-average filtration and yields 
the exponential probability functions given by Eq. Q. 
Therefore we use Eq. ((3|) for estimations of the confi- 
dence probability of peaks in the power spectrum. 

In result we obtain a strong peak of P^(k) at k = 
11 on the significance level well exceeding 5<r. This 
peak manifests the periodicity of the correlation func- 
tion with the period Arj = 0.036 ± 0.002 which is in 
agreement with the periodicity discussed in Sect. [21 



4 Periodicity in distribution of Mg II systems 

It should be noted in this context that a set of sta- 
tistical investigations of resonance absorption doublets 
Mg II (A 2796, 2803 A) on the base of SDSS catalog 
have been performed recently by a few groups (e.g. 
Bouche et aDl2006[|Prochter et al1l2006l . JMenard et al. 
20081 iLundgren et alj|2009h . It follows from their re- 
sults that the spatial distribution of strong doublet 



lines Mg II with W£ 



A2796 



> 



1 A differs from the dis- 



tribution of relatively weak lines W^ 2796 < 1 A, where 
W^ 2796 is the rest frame equivalent width of the line 
A = 2796 A. It seems to be plausible that the strong 
and weak Mg II absorbers are associated with differ- 
ent gala xy populations which display different evolv- 
ing (e.g.jLundgren et al . 2009) and morphological (e.g., 
Kacprzak et al.l 120071 : 120081 1 properties. Therefore one 
could anticipate that radial distributions of both types 
of absorption systems correspond to different tenden- 
cies of clustering and different periodicities. For in- 
stance, weaker absorbers appear to be spreaded over 
space-time more homogeneously and with weaker ef- 
fects of galaxies evolution. 
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Fig. 6 Upper panel: Power spectrum P(k) (5 < k < 50) calcu- 
lated in the point-like statistical approach (Eq. J2}) w ith selection 
of 961 absorption systems Mg II from the catalog bv lBouche et al.l 
(2006) within the interval 0.38 < z < 0.8 (0.35 < n < 0.67) sat- 
isfying a condition of W r < 1.4 A. Lower panel: Comparison of 
two histograms displaying d istributions of 478 Mg I I absorption 
systems from the catalog bv lRvabinkov et al] 1120031) (thick solid 
lines) and 1806 Mg II systems from the catalog bv lBouche et al.l 
(2006) (long dashed lines) with respect to the rest-frame equiva- 
lent width W r . 



Let us note that the main result of the present paper 
is bound up with relatively weak ALSs. It is illustrated 
on the lower panel of Fig. [6] by the histograms display- 
ing the distribution of 478 Mg II absorption doublets 
(A2796, 2803) involved in our analysis with respect to 
W^ 2796 and the analogous distribution o f 1806 Mg II 
systems from the catalog bv lBouche et al.l ( 20061 ) based 
on the third data release (DR3) of the SDSS. The dis- 
tribution of Mg II systems from our catalog has a max- 
imum at w^ 2796 < 0.4 A while the distribution of 



the M g II absorbers from the catalog of iBouche et al 



(J2006I) shows a maximum at W^ 2796 > 1.0 A, i.e., the 
former one is noticeably shifted towards lower values 
of W^ 2796 . Let us note that the more extended cata- 



logue of Mg II systems by iProchter et al.l (|2006l ) also 
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Fig. 7 Upper panel: Power spectrum P(k) calculated with 
using Eq. Q) for 730 QSOs (0.035 < z e < 5.01) from the cata- 
logue |R^abinkgv_eLaL| 1120031 ). whose spectra include 2322 ALSs 
used in the present analysis. An interval of harmonic num- 
bers 25 < k < 70 corresponds to comoving scales (233.3 — 
83.3) /i -1 Mpc, at Q m = 0.23. Lower panel: Same as the up- 
per pa nel but for 96 1 QSO s (0.42 < z e < 2.90) from the cata- 
logue IBouche et al.l 1120061 ), whose spectra include 961 systems 
of Mg II used for the calculations of P(k) in Fig.[S] An interval 
of harmonic numbers 5 < k < 50 corresponds to comoving scales 
(712.9-71.3) h' 1 Mpc. 



based on SDSS data accumulates still stronger lines 
with W^ 2796 > 1 A which are still less consistent with 
our sample of ALSs. As discussed above it is likely that 
these two samples correspond predominantly to differ- 
ent populations of galaxies. 

Nevertheless, we estimate here the presence of the 
same periodicity in the redistribution of Mg II sys- 
tems as discussed in Sections [2] and [3] restricting our- 
selves by relatively weak systems (W^ 2796 < 1 A). The 
upper panel in Fig. |6] represents the power spectrum 
P(/c) calculated according to the point-like approach of 
Eq. E)). We employ N tn t= 961 systems from the cata- 



can see the peak P(fc) at k — 9 for the same period 
A?? = (r/max - Vmin)/k = 0.036 ± 0.002 as in Figs. [Ml 
Additionally, we test possible selection effects which 
might induce periodicities in the initial samples of 
QSOs. Exampl es of such simu lated periodicities are 
given, e.g., by iHartnetd (|2009i ) and iBell fe; Comeau 
(|2010f ). The upper panel in Fig. [7] represents the power 
spectrum P(k) calculated according to Eq. © for 730 
original QSOs within an interval of r/ e — 0.035 — 1.98 
(z e — 0.035 — 5.01) sampled for the statistical treat- 
ment of 2322 ALSs. We see that in the chosen range of 
harmonic numbers, k=25 - 70, there are no significant 
peaks in the power spectrum of the QSOs. Note that 
the period A77 = 0.036 indicated in Sections [2] and [3] 
corresponds to the harmonic number k=54. The lower 
panel demonstrates similar power spectrum calculated 
for 961 QSOs within an interval of r\ e = 0.39 — 1.58 
(z e = 0.42 — 2.90) from the catalog by IBouche et al. 



log of lBoucheetail (J2006I) at W* 2796 < 1.4 A within 
the interval 0.35 < 77 < 0.67 (0.38 < z < 0.80). One 



( 20061 ) . We use their spectra for the test of periodicity 
of 961 relatively weak (W r < 1.4 A ) systems Mg II. 
The power spectrum also does not contain significant 
peaks. Note that the period A77 = 0.036 corresponds to 
the harmonic number k=33. On the other hand, pos- 
sible simulated periodicities at Az e = 0.258, .312, and 
0.44 i n z e -distribution of QSOs indicated by lHartnett 
( 20091 ) for SDSS data correspond to harmonic numbers 
k« 10, 8, and 6, respectively. Additional analysis of 
possible selection effects was carried out in Paper I. 



5 Conclusions and discussion 

On the base of the updated c atalog of absorption 
systems (JRvabinkov et alj 120031) we have performed 
the statistical analysis of 2322 absorption-line systems 
(ALSs) detected in the QSO spectra in the redshift 
range z = 0.0 — 4.3. The ALS redshifts are averaged 
within the velocity interval Sv = 500 km s" 1 along 
each line of sight. This velocity interval containing 
one averaged point corresponds to the spatial scales 
^5 — 6 h^ 1 Mpc (at z = 0.0-4.3) of comoving dis- 
tances for ACDM model at il m — 0.23. The averaged 
points can trace possible large-scale variations of ab- 
sorbing matter. 

The main results can be summarized as follows: 
(1) The distribution of the ALSs relative to the dimen- 
sionless comoving distances r]{z) given by Eq. (jTJ dis- 
plays a statistically significant pattern of alternating 
maxima (peaks) and minima (dips) against a smoother 
curve (trend). This sequence of peaks and dips com- 
prises a certain periodical component with a period of 
Arj = 0.036 ± 0.002. The power spectrum calculated 
according to Eq. ^ for the whole interval of ry (77 = 0.0 



- 1.87) contains the peak at the significance level ex- 
ceeding 3cr relatively to the hypothesis of the Poisson 
^-distribution. In the most statistically representative 
interval rj = 1.148 - 1.870 (z = 1.64 - 4.3) the same pe- 
riodicity reveals itself at the significance level 4.5<r. Still 
more prominent peak corresponding to the same period 
A?y arises in the power spectrum calculated for the two- 
point correlation function £(6r)) with use of Eqs. (0- 
111). 

(2) The same period has been revealed also in the 
power spectra of 961 Mg II (A2796, 2803) abs orption 
systems from the catalog bv lBouche et al.l ( 20061) within 
the interval 0.35 < 77 < 0.67 (0.38 < z < 0.80) calcu- 
lated for relatively small rest frame equivalent widths 
W^ 2796 < 1.4 A. Although the appropriate peak in the 
power spectra is not very significant (slightly exceeds 
3er) this is comparable with the peak obtained in Sec- 
tion [2] for the sample of ALSs within the whole interval 
rj = 0.0 - 1.87 (z = 0.0 - 4.3). It may evidence for avail- 
ability at least a weak periodical component at z < 1 . 

(3) The special feature of the results is the appear- 
ance of the long-range (or intermediate-range) order vi- 
sually displayed by the two-point correlation function 
in Fig. [5] This order corresponds to some stability of 
the phase of periodical component which leads to an in- 
variance of the correlation function under translations 
by a distance multiple of the period A77. The significant 
peak at the lower panel in Fig. [5] may be treated as a 
measure of such an order. 

(4) The dimensionless period A77 obtained here cor- 
responds to the spatial scale of line-of-sight comoving 
distance D c = c/H x A77 = 108 ±6 h^ 1 Mpc. Accord- 
ing to Table Q] this scale D c is consistent with a charac- 
teristic scale of the Large Scale Structure (LSS) in the 
spatial distribution of galaxies and clusters of galaxies 
at relatively small redshifts z < 0.5. On the other hand, 
the appropriate temporal interval T c = 1/H x A?7 = 
350 ± 20 hr 1 Mvr is consi stent with recent results of 
Aref'eva fc Koshelevi lJ2008h . 



Summarized results admit both interpretations of 
the periodicity: spatial one - appearance of partly or- 
dered spatial structures of matter in the early Universe 
(see below) or temporal one - generation of some tem- 
poral wave process in th e course of the cqsm ological 
evolution. For instance, iKaminker et al.l (|2000| ) and 
Paper I argued rather in favour of the temporal inter- 
pretation of the features of ALS distribution. In this 
paper, however, we follow a spatial consideration. Note 
that both hyp otheses have t o obey to the cosmological 
principle (e.g. |Peebleslll993[ ). i.e., partly ordered struc- 
tures or periodical temporal process could be observed 
in any spatial points of the Universe. 

A special wave of interest to the effects of period- 
icity in the redshift distribution of galaxies was initi- 



ated by iBroadhurst et alj ( 1990J). Their pencil-beam 
surveys near the Galactic poles displayed a periodicity 
on a scale about 130 hr 1 Mpc which might be inter- 
preted either as a pure spatial quasi-periodic pattern 

of LS S constituents (set of clumps or walls and voids 
177 — , . „ . — , , I L I r^ 



Kurki-Suonio ct al. 1990, Kaiser & Peacock 



1991 



van de Wevgaertlll991. Dckel ct al. 1992, Yoshida et al 



20011 and references therein) or as a spatial-temporal 
sequence of pronounced and depressed epochs which 
could become apparent in the distribut ion of matter 
(e.g., lMorikawalll99lt Irlirano et alll2008l ). 

Table Q] represents the examples of characteristic 
scales obtained in last years by a few groups of au- 
thors employing statistical analysis of galaxies and 
clusters of galaxies. The characteristic scales in the 
last but one column may be approximately subdivided 
into two groups: the scales belonging to the inter- 
val 70 ± 20 h^ 1 Mpc and the scales within 120 ± 
20 h- 1 Mpc. The period D c = 108 ± 6 /i" 1 Mpc 
gets into the latter group. Esp ecially good co n sent oc - 
curs with results obtain ed b v IPercival et al.l (I2007ah . 
Gaztaiiaga et al.l ( 2009aJ) and iKazin et al.l (|2010al ) for 
the characteristic scale of the Baryon Acoustic Oscilla- 
tions (B AOs) in the large-sca le distribution of matter 
(see also iMartinez et al.ll2009l) . 

Actually, the BAOs are discussed in literature 
in connection with a scale of the sound horizon at 



the epoch of recombinatio n (e.g., iBlake fc Glazebrook 
20031 IPercival et al.ll2007b[) which may be measured as 
a single peculi arity (bump) in the sp atial correlation 
function (e.g., lEisenstein et al.l 120051 ) . On the other 



hand, this spatial scale manifests itself as a series of 
regular variations imprinted on the power spectrum of 
the distribution of galaxies (e.g.. lEisenstein fc Hulll998 , 
Eisenstein et"aHll998l IPercival et al.l l2007bl) . By con- 
trast, we discuss here the regular spatial variations pro- 
viding one prominent peak in the power spectra. In this 
context let us note that the correlation function calcu- 



lated bv IMartinez et al.l (|2009l ) shows a hint of a possi- 
ble secondary peak at a scale of ~ 170 h~ x Mpc. Being- 
confirmed such a tendency would be in consent with 
the hypothesis that at least an intermediate-range or- 
der may be present in the spatial distribution of ALSs. 
Still more pronounced regular spatial variations of 
the standard (see Section[3]) two-point correlation func- 
tion of galax y super cluste rs were obtained, e.g., by 
Einasto et all |l997al) and iTago et al.l (|2002l ). These 
authors indicated the possibility of a long-range (or 
intermediate-range) order with a characteristic period 
of 115 ± 7 h^ 1 Mpc in the spatial distribution of su- 
perclusters: the largest relatively isolated bounded sys- 
tems. Their results are qualitatively consistent with 
the two-point correlation function shown in Fig. [5j 
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Table 1 Characteristic scales of LSS relatively to comoving system in ACDM-model 



Objects 



Surveys 
Catalogs 



Redshifts ACDM Scales References 

fl m h,- 1 Mpc 



92,952 G V SDSS DR1 5 ) (radial) z < 0.14 



0.3 



16,756 G 



LCRS 6 > 



z<0.11 



74 ±17 
60 ±10 



0.15 < z < 0.30 



110.3 ±3.9 



75,000 LRG 



SDSS DR6 



0.25 



14) 



1324 CG 2 ) 


Abell CG 










+ 


+ 


z < 0.28 


ii) 


115 ±7 


15) 


284 X-ray CG 3 > 


RBS 7 ) 










229,193 G 


2dF GRS 8 > 






73.4 ±5.8 




+ 


+ 


z <0.3 


12) 


+ 


16) 


427,512 G 


SDSS DR5 






127 ±21 




143,368 G 


2dF GRS 


z = 0.2 




112.9 ±3.3 




+ 


+ 




0.23 13 > 




17) 


77,801 LRG 4 > 


SDSS DR5 


z = 0.35 




104.7 ±3.2 





18) 







0.40 < z < 0.47 




108.9 ±4.0 




105,831 LRG 


SDSS DR7 


0.16 < z < 0.36 


0.25 


101.7 ±3.0 


19) 


2167 ALS 
961 Mg II 


RKV 9 > 
Bouche 10 > 


< z < 4.3 
0.38 < z < 0.80 


0.23 
0.23 


108 ±6 
108 ±6 


present 
work 



*' G - Galaxies; 2 ' CG - clusters of galaxies; 3 ' X-ray CG - X-ray clusters of galaxies; 4 > LRG - Luminous red galaxies; 5 ' SDSS 
- Sloan Digital Sky Survey; 6 ) LC RS - Las Campanas Red shift Survey; 7 > RBS - R OSAT Bright Survey; 8 ) 2dF GRS - 2dF 
Galaxy Redshift Su rvey; 9 > RKV - iRvabinkov et alj J2003I '); 10) iBouche etldl ||2006| ) ; n) distances betw een clusters were 
calculated using the Mattig] dl958T) formula with the deceleration parameter qo = 0.5 I IEinasto et al.ll200l[) ; 12 ' di stance intervals 
were c alculat ed using Ar = c_Az/Hq; 13 ' the scales r s have bee n recal culated for this Table at f2 m = 0.23; 14 > iDo roshkcvic h et alj 
||2004|'); 15 > iTaeo et alj <2002h : 16 > lHartnett fc Hiranol J2008T ); 17 ) iPercival et al] d2007al 1; 18 ) iGaztanaea et alj J2009arK "*> 
iKazin et alj Il2010al ). Sign + means that the data of both catalogs and both types of objects were used, as well as both scales were 
revealed for the same sample of data. 



Einasto et alj (|l997bl l considered a model of a quasi- 
regular lattice with random distribution of clusters and 
superclusters along lattice edges. They showed that the 
spatial correlation function for a partly regular struc- 
ture displays oscillations similar to those revealed in the 
correlation function of supercluste rs. 

By analogy with Einasto et al.l ( 1997a ) we consider 
in the Appendix two toy models of partly regular struc- 
tures with essential random components and calculate 
the power spectrum of radial distribution of points. The 
power spectra for both models display the peaks with 
approximately the same significance as in Fig. 0J It 
illustrates that rather moderate spatial ordering might 
yield the periodicity of radial distribution similar to the 
sample of ALSs discussed here. 

It should be emphasized that we represent only vi- 
sual examples of partly ordered structure without strict 
arguments in favour of the spatial interpretation of the 



periodicity obtained. This interpretation should be ver- 
ified by special investigations in a future. 
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A Toy models for partly ordered SC lattice 

Figs. |5] and [§] represent two models of simulations 
performed for partly ordered structure of points based 
on the simple-cubic (SC) lattice and the results of ap- 
propriate calculations of power spectra. 

Fig. [5] displays a cloud-like distribution of points 
around the vertices of the SC lattice with unit lat- 
tice constant. We consider a superposition of 3 similar 
lattices with coordinates of the vertices represented as 
x v = i v + Sx v , y v = j v + Sy v , z v — k v + 6z v , where 
i v , j v , and k v are integer numbers within the interval 
from —50 to +50 with respect to the zero point and 
8x v , 8y v , and 5z v are random numbers with the Gaus- 
sian distribution and standard deviation o~ v — 0.3. In 
other words we have a cloud of three random vertices 
instead of one lattice vertex. 

The radial distribution function for the random 
points with radial coordinate r v = \Jx^ + y% + z% 
within a sphere restricted by the radius R — 50 can 
be calculated as a value 



fj, = /j, R (ri) 



Ny(ri) — n inrfAri 
y/no^irrfAri 



(Al) 



where r\ is a mean radius of l—th. concentric layer with 
fixed radial width (bin) An = 0.1, Ny(ri) is a number 
of vertices within a set of non-overlapping layers, no is 
a mean density (number per unit volume) of vertices. 

We calculate the power spectra for each realization 
of [i given by (|A1|) using a version of Eq. ([7]). In this 
version r m substitutes for 6rj m and the radius R stands 
for Sfj ; /i m = /XR,(r m ), where the value of r m runs 
over a uniform set of points (with a step 0.10) r m = 
0.05, 0.15, 0.25, ...49.95 numerated by m = 1,2,...JV P 
and Nfj, = 500. By contrast with Eq. (0 we apply here 
the approach of independent bins (bin=0.1) and put 
U(k) — 1 in the denominator. 

The upper panel in Fig. [5] displays the mean value 
of the power spectra calculated for 100 realizations of 
the radial distribution function (jAlj) . The peak with 
significance ~ 4.5er (cf Fig. [4]) indicates the presence 
of oscillations with a period equal to unity (lattice con- 
stant). 

Similar calculations are represented in Fig. [S] where 
we use so named hole-like lattices constructed on the 
base of SC lattice x v = i v , y v = j v , z v = k v , where 
*-«j jv, and k v are introduced above with a random re- 
moving of some vertices. Parameter u = 0.2 is a ratio of 
the numbers of residual to initial lattice vertices. Addi- 
tionally we immerse the hole-like lattice into the field of 
random Poisson points with sixfold excess of the points 
over the number of residual vertices (e = 6). 

Using the version of Eq. ((7} discussed above, we cal- 
culate the power spectra for different realizations of 



the combined 3-dimensional distribution of points. The 
mean value of the power spectra calculated for 100 real- 
izations of the radial distribution function (|A1|) is rep- 
resented in the upper panel of Fig. [SJ One can see the 
peak with significance ~ 4.5er similar to the case of 
cloud-like lattice. 

The lower panels in Figs. [5] and H] represent 2- 
dimensional illustrations of the partly ordered struc- 
tures used in the power-spectrum calculations. Let us 
emphasize a low degree of visual ordering of the struc- 
tures discussed. 
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Fig. 8 Upper panel: Mean value of 100 realizations of power 
spectra calculated for the radial distribution function Eq. JAll l 
in the case of a cloud-like lattice (see text) constructed as su- 
perposition of 3 random cubic lattices with the standard devi- 
ation tTp = 0.3 from simple cubic. Lower panel: fragment of 
2-dimensional illustration displayed one of 100 realizations of the 
cloud-like lattice (calculated for the upper panel); points repre- 
sent lattice vertices and the filled circle indicates the center of 
the sphere. 
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Fig. 9 Same as in Fig. [8] but for superposition of a hole- 
like lattice with a parameter u = 0.2 (see text) and the field of 
random Poisson points with sixfold (e = 6) excess of random 
points over hole-like lattice vertices. Upper panel: mean value of 
100 realization of power spectra; lower panel: filled circles stand 
for the vertices of the hole-like lattice. 



